Spatial quantum correlations in multiple scattered light 



in 
o 
o 

> 

o 
o 



m ■ 

> : 

m . 
m ■ 
O ■ 

<n : 
o . 
vn ■ 
o ■ 

^ : 

a: 
i 

^— > ■ 

G ■ 
cd 

3 : 



P. LodahlQ A.P. Mosk, and A. Lagendijl|f| 

Complex Photonic Systems, Department of Science and Technology and MESA + research institute, 
University of Twente, P.O. Box 217, 7500 AE Enschede, The Netherlands. 
(Dated: February 1, 2008) 

We predict a new spatial quantum correlation in light propagating through a multiple scattering 
random medium. The correlation depends on the quantum state of the light illuminating the 
medium, is infinite range, and dominates over classical mesoscopic intensity correlations. The spatial 
quantum correlation is revealed in the quantum fluctuations of the total transmission or reflection 
through the sample and should be readily observable experimentally. 

PACS numbers: 42.25.Dd, 42.50.Lc, 78.67.-n 



When light propagates through a disordered scatter- 
ing medium, a strongly modulated interference structure, 
known as a speckle pattern, is generated. Speckle pat- 
terns are the most direct manifestation of wave coherence 
in transport of light through samples that are thicker 
than the transport mean free path £, which is the aver- 
age distance over which the direction of light is diffused 
due to random scattering. By analyzing the statistical 
properties, such volume speckle patterns reveal strong 
correlations that are responsible for fundamental physi- 
cal phenomena as the memory effect yj, and enhanced 
mesoscopic fluctuations |3j, Lit LJ III LiJ • Furthermore, clear 
signatures of Anderson localization of light have been ob- 
served by analyzing intensity fluctuations §. 

While spatial correlations of the intensity of multiply 
scattered light have been investigated intensely, spatial 
correlations hidden in the quantum fluctuations of light 
have to our knowledge not been addressed. Here we pre- 
dict that strong quantum correlations exist between spa- 
tially separated parts of a far-field speckle pattern that 
have no classical analog. The quantum correlation is in- 
finite range, since it is independent of the angular sepa- 
ration between the parts of the speckle pattern under in- 
vestigation, and dominates over mesoscopic fluctuations. 
Therefore speckle patterns generated with quantum light 
are much stronger correlated than classical theory pre- 
dicts, and the quantum corrections should be observable 
even for moderate amount of scattering. In contrast, 
the magnitude of classical intensity correlations are all 
found to scale with the inverse of the mesoscopic conduc- 
tance g = N£/L where N is the number of conducting 
modes and L is the thickness of the scattering medium 
0,0,0, 0,0|- Only for extremely strong scattering, i.e. 
close to the Anderson localization transition (g ~ 1), do 
such mesoscopic correlations dominate. 

The quantum correlation will be shown to depend on 
the quantum state of light illuminating the multiple scat- 
tering sample. The physical origin of the effect is due 
to distribution of photon correlations of the input state 
over the spatial degrees of freedom of the speckle pattern. 
This mechanism was previously predicted to lead to new 
long-range correlations of thermal radiation from a dis- 



ordered waveguide Jjj . We evaluate the correlation func- 
tion for three different quantum states: coherent state, 
thermal state, and Fock state. For a coherent state the 
spatial correlation vanishes, while strong correlations and 
anti-correlations are found for the thermal state and Fock 
state, respectively. We show that the spatial quantum 
correlation can be measured conveniently by recording 
the quantum noise of the total transmission or reflection 
from the medium. Such measurements were carried out 
very recently for a coherent state |lfj|. and appear feasi- 
ble also in the near future for other quantum states. Our 
work is connected to the field of quantum imaging where 
spatial correlations are employed for reconstructing im- 
ages 

We describe the propagation of quantum states of light 
through a non-absorbing, multiple scattering medium us- 
ing the formalism developed by Beenakker et al. \vX\\^ . 
The model describes effectively a quasi-lD configuration, 
but this approximation is known to be excellent also for 
a 3D slab geometry. The annihilation operator a b for an 
output mode b is coupled to the annihilation operators 
associated with all input modes and a™ through the 
relation a b — J2 a ' Vfc^o™ + Y^b> r b'ba b ?, where t and r are 
electric field transmission and reflection coefficients, re- 
spectively, and the summations are over all N input and 
output modes. We use the notation that a denotes modes 
to the left and b modes to the right of the medium, c.f. 
Fig. n The modes correspond to different propagation 
directions (k vectors) of the incoming and outgoing light 
as measured directly in the far-field. 

In the following, we explicitly outline the calculation 
of the quantum fluctuations of the light transmitted 
through the sample when light is coupled through a sin- 
gle input mode a. In the quantum description vacuum 
fluctuations coupled into all other input modes a' a 
and b' must be included. The operator n b = a\a b de- 
scribes the number of photons transmitted to the output 
mode b, and the total number of transmitted photons is 
obtained by adding all output modes: ut = J2 b fib- In an 
experiment, the total transmission can be measured with 
an integrating sphere, cf. Fig. ^ The quantum fluctu- 
ations arc quantified through the variance of the photon 
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FIG. 1: (color online). Total transmission through a multi- 
ple scattering sample. A quantum state of light is coupled 
through mode a and the total number of photons is collected 
with an integrating sphere. Vacuum fluctuations in all open 
channels a 1 7^ a and b' must be included in the quantum de- 
scription. 



number: An 2 = (n 2 ) — (n) 2 , where (•) is the quantum 
mechanical expectation value. The fluctuations of the 
total transmission are straightforwardly found to be 

A? 4 = E An fc + E E ( h b n bl ) - (n bo ) (n bl ) , (1) 

b bo 61 /bo 

where an index specifying the fixed input mode (a) has 
been omitted for brevity of notation. We immediately 
see that the variance of the total transmission is differ- 
ent from the the sum of the variances of all outgoing 
modes. Also cross-correlations between different output 
modes contribute, and these correlations are shown to be 
significant in the following. 

In order to evaluate Eq. (JTJ, (n b ) and (h bo h bl ) must 
be computed. The non-vanishing contributions are 
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(2a) 
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The expectation values of the various operator products can be evaluated using that vacuum is coupled through all 
other channels than a. We obtain 



a a' a b' a b" a a" 



- (flT) &a> ,aS a " ,afo' ,b" 



(3a) 
(3b) 
(3c) 



By combining Eqs. |J2J and using an identity de- 
rived from commutation relations the photon num- 
ber variance of mode b and the cross-correlation between 
two different directions bo ^ b\ can be found 



= T ab + T a b 2 (F™ - 1) (2 + 8/3. g ) , (4a) 



(h bo h bl ) - (h b0 ) (rW = —2 {F , n + 4/3ff) 



(4b) 



Here we have introduced the intensity transmission co- 
efficients T a b — \tab\ 2 that have been averaged over all 
realizations of disorder (denoted by a bar). Products of 
transmission coefficients are expanded to first-order in 
l/g, i.e. T abo T abl « (1 + S baM ) (1 + 4/3 ff )7^ 2 15]. The 



Fano factor F™ — A(n™) 2 / (n™) measures the variance 
of the photon number relative to the average number 
of photons, and is equal to unity for Poissonian pho- 
ton statistics. We evaluate the Fano factor for three 
different single- mode quantum states of light: a coher- 
ent state (CS), a thermal state (TS), and an n-photon 
Fock state (FS), which correspond to (F™) cs = 1, 
(F l a n ) TS = 1 + (nf) and (F™) FS = 0, respectively 0. 

The classical transmission coefficient can be calculated 
from the theory of light diffusion, which leads to T a b — 
£/NL in the absence of absorption ^t|. Inserting Eqs. 
Q in Eq. QJ, we arrive at the final expression for the 
quantum fluctuations of the total transmission 
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Neglecting at present correction terms of order 1/g, the 
variance of the total transmission is seen to contain a 
linear and a quadratic term in i/L. The former is the 
sum of the variances of each individual channel, which 
corresponds to incoherently adding the fluctuations of 
each output mode. The quadratic term is due to correla- 
tions between different output modes. Note that while a 
quadratic scaling was found also for transmission of clas- 
sical noise through a multiple scattering medium 
the contribution discussed here is a quantum effect. The 
fluctuations of the total reflection can be calculated in a 
similar way, and neglecting correction terms of order 1 jg 
leads to 
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Figure displays the variance of the total transmis- 
sion and reflection for three different quantum states of 
light as a function of the ratio of the mean free path i to 
the sample length L. For a coherent state, the variance 
scales linearly, which was experimentally confirmed in 
|lfj|. For other quantum states of light, a quadratic term 
contributes by either a positive (for the thermal state) or 
negative amount (for the Fock state), which means that 
the photon statistics of the total transmission and reflec- 
tion is super-Poissonian or sub-Poissonian, respectively. 
These quantum correlations are found to be most pro- 
nounced in the reflection for thick samples (i/L — » 0) and 
in the transmission for thin samples (i/L — > I ). This be- 
havior can be understood intuitively since for thick (thin) 
samples relatively few photons are transmitted (reflected) 
and vacuum fluctuations coupled through all open modes 
lead to Poissonian photon statistics. The validity of the 
diffusion model for light propagation requires i/L < i, 
and in this limit each transmitted or reflected photon 
has experienced numerous scattering events. Remark- 
ably the quantum correlations of these multiply scattered 
photons are preserved, and for i/L — > 0, the fluctuations 
of the total reflection equal the fluctuations of the input 
light. This clearly demonstrates that deep in the multi- 
ple scattering regime, the quantum corrections are most 
pronounced in the reflection, and thereby devises a route 
for experimental observation of quantum correlations in 
a volume speckle pattern. 

The quadratic terms in Eqs. (J5J and JfJJl originate from 
spatial correlations between different output directions 
in the speckle pattern. We define the two-point spatial 
correlation function, which is evaluated using Eqs. J3J), 
© and 



(h bo ) (n bl ) 



= 1 



(7) 



where bo ^ b\. Cb b 1 gauges the correlation between pho- 
tons in modes bo and b\ [f8| and depends only on the 
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FIG. 2: (color online). Variance of the quantum fluctuations 
of the total reflection (upper plot) and total transmission 
(lower plot) for a coherent state (full line), a thermal state 
with (n™) = f (dash-dotted line) and an arbitrary Fock state 
(dotted line) as a function of l/L. All correction terms of or- 
der 1/g have been neglected. In the reflection (transmission) 
pronounced variations in the photon statistics are found for 
decreasing (increasing) i/L. 



statistics of the input light. Note that Eq. is ex- 
act and independent of 1/g. The correlation function dis- 
cussed here is fundamentally different from the one en- 
countered in simpler systems as, e.g., a dualport beam- 
splitter, since Eq. Q is averaged over all ensembles of 
disorder. To measure the two-point correlation function 
generally requires involved coincidence detection schemes 
between different output directions from the medium. 
However, we find here that the quantum fluctuations of 
the total transmission, or equivalently the total reflec- 
tion, provide an alternative way of extracting the two- 
point correlation. The total transmission and the quan- 
tum correlation function are connected through 
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where again corrections of order 1/g have been neglected. 

Figure |21 shows the two-point correlation function for 
three different quantum states when varying the average 
number of photons in the input state. For a coherent 
state with an arbitrary average number of photons the 
correlation function is unity, i.e. the spatial correlations 
vanish. For single mode thermal states, the spatial cor- 
relation function is always equal to 2, which means that 
different parts of the speckle pattern are strongly corre- 
lated. For Fock states, a striking non-classical behavior is 
found: the spatial correlation function is reduced below 
unity, which means that the spatial directions are anti- 
correlated. This corresponds to spatial anti-bunching 
of light, and the correlation function Cb obl is the spa- 
tial counterpart of the second-order coherence function 
g^(r) that is omnipresent in quantum optics [fij . For 
a single-photon Fock state, the correlation function van- 
ishes identically since if a photon is detected in one chan- 
nel the probability of detecting another photon in a dif- 
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FIG. 3: (color online). Two-point correlation function Cb bi 
for three different quantum states of light as a function of the 
average number of photons in the incoming field. The cor- 
relation function exhibits either strong correlation (thermal 
state) or anti-correlation (Fock state) relative to the uncorre- 
cted case (coherent state). 



ferent channel is zero. For Fock states with an increasing 
number of photons, the correlation function approaches 
unity, i.e. the spatial correlations vanish. 

In addition to the spatial correlations also mesoscopic 
fluctuations contribute to the quantum fluctuations of 
the total transmission, c.f. Eq. JSJ). Figure 0] displays 
these fluctuations as a function of the mesoscopic con- 
ductance g. For small values of g, the mesoscopic corre- 
lations either enhance or reduce the quantum fluctuations 
for thermal states and Fock states, respectively. In con- 
trast, the quantum fluctuations for a coherent state of 
light are unaltered. These mesoscopic fluctuations pro- 
vide another example of correlations that are present only 
in a quantum description of multiple scattering. 

It is instructive to compare our findings on photon fluc- 
tuations in diffusive propagation to the transport of elec- 
trons through an elastically scattering disordered conduc- 
tor. In the diffusive regime of electron transport, Poisso- 
nian electron fluctuations are reduced, which is a direct 
consequence of the Fermi-Dirac statistics of electrons in 
each individual conduction channel [20ll2l| . However, no 
correlations exist between different conduction channels 
and the fluctuations can simply be added up incoherently. 
Our work demonstrates that the dominating correlation 
in diffusive transport of photons is due to correlation be- 
tween different spatial modes that are the optical anal- 
ogy to electronic conduction channels, thus pointing out 
a fundamental difference between electrons and photons. 

We have predicted a spatial correlation present in a 
quantum description of multiple scattering of light. The 
spatial correlation depends on the quantum state of the 
light source illuminating the sample, and can be mea- 
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FIG. 4: (color online). Fluctuations of the total transmission 
for a coherent state, a thermal state with (n^T) = 1, and an 
arbitrary Fock state as a function of the mesoscopic conduc- 
tance g, and for fixed amount of disorder: l/L = 1/3. For 
decreasing g, the mesoscopic correlations either enhance or 
suppress quantum fluctuations for the thermal state and Fock 
state, respectively. The fluctuations of the coherent state are 
unaltered for all values of g. When g approaches unity, the 
localization regime is entered and the model is invalid. 



sured by recording the quantum fluctuations of the total 
transmission or reflection from the sample. Our work 
shows that multiple scattering of non-classical light in- 
duces stronger correlations than for classical light. The 
predictions should be readily observable experimentally 
even for moderate amounts of scattering. 
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"Stichting voor Fundamentcel Onderzoek der Materie 
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